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We calculate the second order roughness correction to the Casimir energy for two parallel metallic 
mirrors. Our results may also be applied to the plane-sphere geometry used in most experiments. 
The metallic mirrors are described by the plasma model, with arbitrary values for the plasma wave- 
length, the mirror separation and the roughness correlation length, with the roughness amplitude 
remaining the smallest length scale for perturbation theory to hold. From the analysis of the intra- 
cavity field fluctuations, we obtain the Casimir energy correction in terms of generalized reflection 
operators, which account for diffraction and polarization coupling in the scattering by the rough 
surfaces. We present simple analytical expressions for several limiting cases, as well as numerical 
results that allow for a reliable calculation of the roughness correction in real experiments. The 
correction is larger than the result of the Proximity Force Approximation, which is obtained from 
our theory as a limiting case (very smooth surfaces). 

PACS numbers: 42.50.-p, 03.70.+k, 68.35.Ct 
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I. INTRODUCTION 

The Casimir force of attraction between metallic mir- 
rors Q] has been measured with high experimental pre- 
cision over the last few years . These new experiments 
allow for an accurate theory/experiment comparison 
opening the way for the search for new weak forces with 
submillimetric ranges Q. On the theoretical front, accu- 
rate results based on realistic models are sorely needed in 
order to match the desired levels of accuracy. Three im- 
portant effects provide the main corrections to the ideal 
configuration considered by Casimir: non zero temper- 
ature 0, finite conductivity 0, Q and roughness of the 
mirrors |8j|-|ll|. Temperature corrections are important 
when the distance L between the mirrors is above 1/zm, 
whereas finite conductivity and roughness provide the 
major corrections for the short distances (of the order of 
a few hundred nanometers) probed by most experiments. 

In principle, these effects must be taken into account 
simultaneously. The overall correction is not in general 
the product of the separate corrections calculated inde- 
pendently. In particular, the correlation between finite 
conductivity and roughness effects is essential, because 
they both intervene at the same range of L. Therefore, 
a reliable theory for short distances must analyze the 
roughness effect in the context of a finite-conductivity 
model for the material medium. To this aim, we de- 
scribe the optical properties of the metallic mirrors by 
the plasma model. 

When the surface profiles are nearly smooth over dis- 
tances of the order of L, the roughness correction may 
be calculated from the Proximity Force Approximation 
(PFA) 12]. In this approximation, the Casimir energy is 
computed from the formula for parallel planes by aver- 
aging the 'local' distance over the surface In order 
to derive more general results, we develop a perturbative 
theory for the Casimir energy with rough plane mirrors, 



allowing for the computation of the energy correction 
when the surface profile varies on arbitrarily short length 
scales, provided that they are larger than the rough- 
ness amplitude (otherwise the perturbative approxima- 
tion would not apply). Our approach is applicable to 
most Casimir force measurements between metallic mir- 
rors. 



We follow the approach of Ref. 0], and consider 
the two mirrors as a plane Fabry-Perot cavity, which is 
treated as composed optical network in order to calcu- 
late the intracavity field fluctuations. We then derive a 
formal result for the Casimir energy up to second order 
in the amplitude of roughness, in terms of generalized 
reflection coefficients describing the scattering by rough 
surfaces, taking into account the coupling between Trans- 
verse Electric (TE) and Transverse Magnetic (TM) po- 
larizations. Numerical results derived from the present 
calculation were presented in a letter [l4|, together with 
some analytical limiting cases. In this paper we present 
the complete derivation and the explicit formulas used in 
0. 



This paper is organized in the following way. In Sec. II, 
we present some basic definitions and assumptions, and 
discuss the validity of the PFA in two different contexts. 
In Sec. Ill, we derive the formal, general result for the sec- 
ond order roughness energy correction, which is then ap- 
plied to the specific plasma-model calculation presented 
in Sec. IV. Several limiting cases are discussed in the fol- 
lowing sections: the short roughness wavelength regime 
(Sec. V), the perfectly- reflecting limit (Sec. VI) and the 
plasmon limit (Sec. VII). In Sec. VIII, we discuss the 
example of a Gaussian roughness spectrum and present 
some concluding remarks. Three appendices present ad- 
ditional details of the derivations. 
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FIG. 1: Fabry-Perot cavity of length L. 



GENERAL CONSIDERATIONS AND 
ASSUMPTIONS 



Our Fabry-Perot cavity of length L is composed of 
two parallel mirrors with rough surfaces, as shown in 
Fig. 1. We analyze the cavity as a composed optical 
network, and calculate the fluctuations of the intracavity 
fields propagating along the positive and negative z-axis, 

Ec an d Eci m terms of the fluctuations of the incom- 
ing free-space fields and EJ£ (also shown in Fig. 1 
are the outgoing fields E£ ut and E™'). In appendix A, 
we show that the Casimir energy turns out to depend 
only on the coefficients describing the reflection of the 
intracavity fields by the internal sides of mirrors Ml and 
M2. The functions hi (r) and /12 (r) define their surface 
profiles with respect to reference planes at z = (see 
Fig. 2) and z — L, respectively, r collects the two trans- 
verse coordinates (x, y) orthogonal to the cavity exten- 
sion. By construction, both h\ and hi have zero spatial 
averages: (hj) = 0,j = 1,2, and are counted as positive 
when they correspond to local length decreases below the 
mean value L. 

We assume that the two surfaces are statistically inde- 
pendent, so that the cross correlation function vanishes: 



(fti(r)ft 2 (r')) = 0. 



(1) 



Translational symmetry on the xy plane implies that the 
self correlation functions satisfy 

{h j {v)h j {v'))={h j {v-v')hM), 3 = E2. 

Then in the Fourier domain we have (k is a two- 
dimensional vector) 



(^■(k)^-(k')) = (2 7 r) 2 <5( 2 '(k + k')a jj (k) 



(2) 



where Hj(k) is the Fourier transformation of hj(r) and 
the roughness spectrum Ujj (k) is the Fourier transform 
of the self correlation function: 

oa(k) = J d 2 re- lk ' r (Mr)M0)). 



h( x >y) 



z = z 
FIG. 2: Magnified detail of the internal surface of mirror Ml. 

We assume that the area A of the mirrors contains 
many correlation areas: A ^> £q, where £q is the correla- 
tion length characteristic of the self correlation function. 
In this case, a single mirror already contains many inde- 
pendent realizations of surface profiles, and hence spatial 
and ensemble averages are equivalent. 

We also assume that the deformation amplitudes are 
very small, in the scale of the mean cavity length, \hj\ <C 
L, as well as in the scale of the correlation length, \hj\ <C 
Ic, so that the surface profile gradients satisfy |V/ij| <C 
1. This allows us to treat the surface deformations as 
small perturbations of the ideal plane geometry. The 
Casimir energy is then calculated up to second order of 
the deformation amplitudes. More precisely, the energy 
correction is obtained in terms of 



(iJ i (k)^(-k))=A^ i (k), 



(3) 



where we have used Eq. (J2J. 

Most experiments are performed with a plane-sphere 
(PS) setup, instead of the plane-plane (PP) cavity used 
as the benchmark for our perturbative calculation. How- 
ever, our results may also be applied to those experi- 
ments, provided that we use the PFA to connect the two 
different geometries. In this case, the force Fpg between 
a sphere of radius R and a plane at a distance of closest 
approach L is given in terms of the energy -Epp for the 
plane-plane cavity as follows: 



Fp S (L) = 2-kR 



Epp (£) 
A 



(4) 



The relative roughness correction of the force in the 
plane-sphere geometry may then be obtained from the 
relative energy correction calculated in this paper: 



SF, 



PS 



5E } 



pp 



Fps Epp 



(5) 



We emphasize that the PFA amounts to the addition of 
contributions corresponding to different local inter-plate 
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distances, assuming these contributions to be indepen- 
dent. But the Casimir energy is not additive, so that 
the PFA cannot be exact, although it is often improperly 
called a theorem. 

Note that the conditions required for applying the PFA 
to connect the plane-sphere and the plane-plane geome- 
tries are quite different from those necessary for using 
the PFA in the computation of the roughness correction 
itself. In the first case, it is necessary that the radius R 
is large enough, so that the separation L satisfies L -C R 
|15|. Moreover, to avoid any interplay between curvature 
and roughness effects, one requires the correlation length 
lc to be small enough, so that many correlation areas 
are contained in a given nearly-plane local section of the 
spherical surface: £q <C RL. In contrast, applying the 
PFA to roughness requires the surfaces to be nearly plane 
in the scale of the separation: £c ^> L. Then, the fol- 
lowing second-order roughness correction to the Casimir 
energy is obtained [TTI | : 

6E PP *^±(h 2 + h 2 ), (PFA) (6) 

Hence, the PFA result for the roughness correction de- 
pends on the second order derivative of the energy and 
on the variances of the length deformations hi and h 2 . 
This expression is equivalent to the procedure used for 
analyzing the effect of roughness in recent experiments 
0, In the following sections, we will assume the 

PFA to provide a valid description of curvature, but not 
of roughness. 

III. REFLECTION AND LOOP FUNCTIONS 
FOR ROUGH MIRRORS 

In this section we start to develop a description of re- 
flection by the internal sides of the cavity mirrors (see 
Fig. 1), leaving the more general theory which takes into 
account the coupling with the external fields to appendix 
A. 

We take the mixed Fourier representation for the in- 
tracavity fields: 

— > _>TE _>TM 

Ec (k, z, u) = (E c (k, W )e TE + E c (k, w)e™) e^ z , 

(7) 

where u> is the frequency, k is the two-dimensional 
wavevector associated to propagation parallel to the xy 
plane, and k z = sgn(w) y/ 'oj 2 / 'c 2 — k 2 , with sgn denoting 
the sign function. The complete wavevector is given by 

K = k + k z z. The field Ec is written in a similar way, 
except for the replacement k z — ► — k z . The TE and TM 
unitary vectors are defined in the following way: 

e TE = 2 x k, (8) 



e™ = e TE x K. (9) 



It is useful to employ the Dirac notation, with the ket 

| Ec 0*0) providing a compact notation for the field am- 
plitudes: 

E c (k,w) = (k,p| Ec 0*0), 

where p = TE, TM denotes the polarization. 

The reflection by the mirror Ml at frequency u> is writ- 
ten as 

| Ec 0*0) = Kx(u)\ Ec 0*0)- (10) 

When the surface is rough, the operator 1Zi(u>) mixes 
up different values of k and polarizations. On the other 
hand, the field frequency is conserved, since the surface 
is at rest. The explicit form of 1|10[1 is 

Eo (k,c*0 = / ^^(M^MIkV) Eo (k»- 

(11) 

The reflection by mirror M2 is defined in a similar way 
in terms of the operator 72-2 0*0' 

We expand the reflection operators TZj(w), j = 1, 2, in 
powers of the deformation amplitudes hj : 

Tij ( w ) = nf («) + snf («) + snf ( w ). (12) 

The zero-th order operators 72-^ correspond to ideally 
plane surfaces. They do not modify the polarization nor 
the momentum k, and hence are diagonal in the basis 
{IM>: 

(k,p|7efM|kV) = (2vr)V 2 )(k-k') V^(k^), 

where rj(k, ui) are the specular reflection coefficients for 
a plane mirror. 

For the ideal Fabry-Perot cavity, the Casimir effect 
may be entirely described by these reflection coefficients, 
which characterize the optical properties of the cavity 
as seen by the intracavity field [17|, As shown in Ap- 
pendix A, a similar result holds for a cavity with rough 
mirrors, except that the specular reflection coefficients 
are replaced by the reflection operators defined above. 
The Casimir force is calculated from the spectral den- 
sity characterizing the vacuum field fluctuations. For the 
intracavity field, the free-space spectral density for polar- 
ization p is multiplied by the generalized Airy function 
g p (k, u>), which quantifies the joint boundary effect of the 
two mirrors. We then derive the Casimir force after in- 
cluding the contribution of evanescent waves: 

Fpp is defined as the z-component of the force on mirror 
Ml; hence it is positive in case of attraction. 

We compute <? p (k, ui) up to second order in hi and hi : 

ffp (k,w) = . 9 (°)(k, w )+ ( 5 5 W(k, w )+ ( 5 5 ( 2 )(k, w ). (15) 
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1 is the Airy function for the ideal plane cav- 



9i> 
ity [1 



4 >(k,w)=l + / p (k,w) + / p (k,u;)* 



where / p (k, w) is the corresponding loop function. It is 
given by the superposition of all propagation factors rep- 
resenting a closed loop with n round-trips inside the cav- 
ity: 



/ P (k,w) 



n=l 



(r?(k jW )rf(k jW ) 



-2kL\ 



r?(k,o;)rf(k,a;)e- 2 ' cL 
l-r?(k,w)rf(k,u;)e- 2K - L: 



(16) 



where «; = —i^Lu 2 /c 2 — k 2 . When replacing g p by g p in 
(|14f> . we find the well-known result for the Casimir force 
in the ideal case [Ti|. 

The first-order Casimir force correction, coming from 
5gp (k,o;) in (|15|1 . vanishes because it is proportional to 
the averages (/ii) and (fo)- Thus, the roughness correc- 
tion is of second order, and results from the contribu- 
tion of Sg p (k, uj). These functions are written in terms 



of 'rough' loop functions Sfp 2 (k, uj) and Sf p 



(2ii), 



k,wl 



gathering the second-order contributions of the first 



(81Zj ) and second order (JTS.V^) reflection operators, re- 
spectively: 

5gW (k, c) = 6$ (k, W ) + Sfj 2 «) (k, «) + c. 



(17) 



(5/p 2 (k, uj) is the superposition of all closed loops in- 
volving a single second-order rough reflection at one of 
the mirrors: 



i piv^y 1 ^ ( W ) e -^ )L <wzf ( W ) PH- 1 |k, P ) , 



(18) 



r 



with [j + 1] representing a sum modulo 2. Like 7^ 0) ( w ) 
in Eq. (|13fl . the operators £>(w) and /C(w) are diagonal 
operators, with elements 1 — ff(k, w)r^ (k, w)e~ 2KL and 
re, respectively. 

To understand why <5/ p 2lI ^(k, w) is a generalization of 
the ideal loop function / p (k, uj), we should read the r.-h.- 
s. of Eq. l|18|l from right to left. The second-order rough 
reflection at mirror j is followed by a one-way propaga- 
tion between the two mirrors (operator exp(— tC(uj)L)), 
and then by a specular reflection at mirror [j + 1]. The 
loop is closed by a second one-way propagation back to 
mirror j. This loop can be preceeded and/or followed 
by arbitrary numbers of round-trips with specular reflec- 
tions, hence the entire expression is sandwiched between 
two operators X>(£j) _1 . 



J 



'Closing the loop' means to ensure that the initial and 
final states are the same, which is represented by the 
ket \k,p) and the corresponding bra in Eq. I|18l) . Since 
all zero-th order processes conserve momentum and po- 
larization, only second-order rough reflections that also 
conserve momentum and polarization are allowed, so 

(2) 

that only diagonal elements of S1Z; are expected to 



contribute in Eq. 
yields 



ljT%|) . Its explicit evaluation indeed 



3 = 1 V 1 



rf. +1] (k,o,)(M«ef M|k,p)e 



(k,tj)rf(k,w)e- 2KL ) 2 



Those diagonal matrix elements are of the form 



(19) 



k,p\6nf\uj)\k,p) = / ^i?g(k,k';c)|^(k-k')| 2 , 



(20) 



where the non-specular coefficients _fc 2 '(k, k'; to) are in- dependent of the profile functions Hj(k). 
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On the other hand, nondiagonal matrix elements of 



5K 



(i) 



contribute to the loop function Sfp ^(k, oj), be- 



cause the latter contains two first-order rough reflections 
instead of just one second-order reflection. These ele- 
ments are of the form 

(MdftfMlk'.p') =iZ^(k 1 k / ; W )fl J -(k-k / ), (21) 

where again the coefficients R^Li (k, k';w) are indepen- 
dent of Hj(k). According to this expression, a given 
Fourier component Ak of the surface profile leads to a 
field momentum modification by Ak. 

Since the elements (k,p|<J7?i- (uj)\k.' , p') are propor- 



tional to Hj(k — k'), the terms associated to first-order 
rough reflections at different mirrors are proportional ei- 
ther to the product 



ffi(k-k')#2(k'-k) 



or to its complex conjugate. As discussed in Sec. II, we 
assume that their average values vanish because the two 
surface profiles are statistically independent. Thus, we 
only keep the terms associated to two first-order rough 
reflections at the same mirror when deriving Sfp 2 ^ (k, w) 
(omitting the dependence with u) in the r.-h.-s.): 



*/f"(k,a>) 



3=1 



k,p\V~ 



-K.L 



TZ (0) e 



l e -KL 



K 



(o) 

[3+1] 



-KL 



8nf ] v- l Wp). 



(22) 



r 



As for Eq. . the sequence of events associated to 
these loops may be read from right to left in the r.-h.-s. 
of (1221) : first-order rough reflection by mirror j , one-way 
propagation to mirror [j + 1], specular reflection, one- 
way back to j, second first-order rough reflection by j. To 



restore the initial sense of propagation, the loop is closed 
by one or more specular round-trips. Before each rough 
reflection, arbitrary numbers of specular round-trips are 
allowed. Explicit evaluation yields 



(2i) 



3=1 V' K ' 



d 2 k' e-^+^ L r? 1 (k)rP; i (k')<k,p|^( 1) (c)|k',p')(k',p'|^( 1) (c)|k,p) 



(1 



rf (k, w)rf (k, uj)e-^ L f (1 - rf (k', w)rf ' (k', uo)e- 2K ' L ) 
I 



(23) 



where k' — —i^/uj 2 /c 2 — k' 2 . 

Before replacing all these results into l|14fl . we first 
write the Casimir force as the real part of integrals of 
the loop functions [l^. Since these functions are ana- 
lytical, Cauchy theorem allows us to replace the integral 
over real frequencies by an integral over the imaginary 
axis in the complex plane of frequency. As a result, lo 
is replaced by £ = — iu, and exp(— kL) becomes a real 



exponential factor, with 



k 2 + J > 0. 



The resulting integrals turn out to be real, and the rough- 
ness correction is then given by 



SFpp(L) = 2 A 



d 2 k 



—Hk 
2tt ^ 

p 



[<54 2i )(k,C) + <54 2ii )(k,0 



(24) 



r 



According to JTHJ and QDJl, <5/i 2ii) (k, £) is given by 
an integral over k' with the integrand proportional to 
|-£/j(k— k')| 2 . We also obtain this factor when computing 



5f^ 2i) (k^) from J23J and (J2HJ. According to ©, when 
averaged it yields A<jjj (k — k'), due to translational sym- 
metry on the xy plane. Hence both loop functions are 
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independent of A, yielding a force proportional to A as 
expected. 

The energy correction is computed from Ij24(l by a sim- 
ple integration: 

/oo 
SF PP (L')dL'. 

In order to simplify the notation, we consider two mirrors 
made of the same metal, and hence with the same optical 
properties (otherwise the correction is given by a trivial 
extension) . Then, after transforming k — k' into k by a 
trivial change of integration variable and taking Q into 
account, we hnd 

SEp P = J ^G(k)a(k), (25) 
I 



with cr(k) = <7ii(k) + er 2 2(k). As discussed in connection 
with Eq. I|21|). the k in G(k) represents the field momen- 
tum transfer induced by a given Fourier component of 
the surface profile. The second-order roughness response 
function G(k) is given by the following general expression 

oo 

G(k) =" M /l/s 6k '' k '- k (26) 



The contribution of the first-order reflection operator is 
calculated from Eq. (|23[1 (from now on we drop the index 
j indicating one of the two mirrors): 



1 e -2(«'+«")L r p' (k / 5 rP " (k // ; ^ R W n ( £. k / ; k // } R W pi ( £. k // j k / 



E 

p'p" 



rP'(k',£)"e 



2 — 2Kj'^j 



) (1 - r*'' 



(k",^) 2 e- 2K " 



(28) 



whereas Eq. i|19|) leads to the following contribution from 
the second-order operator: 



e- 2re ' i r p (k / ,e)i4 ) (C;k',k") 



l-rP(k',£) 2 e 



-2k'L 



(29) 



The response function G(k) is entirely determined by 

the mirrors' non-specular coefficients R^2 and to- 
gether with the specular reflection coefficients and the 
exponential factors describing round-trip propagation in- 
side the cavity. For isotropic material media, symme- 
try requires the response function to depend only on the 
modulus k — |k|. According to this k dependence 

describes the spectral sensitivity of the Casimir energy to 
roughness. Hence, in general the Casimir energy depends 
on the details of the roughness spectrum c(k), and not 
only on the roughness variance 



(hi 



h%) = 



d 2 k 

(27T) 2 



<7(k). 



When G(k) is known and the roughness spectrum mea- 
sured experimentally, Eq. (|25|l allows for a precise and 
straightforward calculation of the roughness correction 
to the Casimir force without using the proximity force 
approximation. 

The PEA is recovered only when the the surface is very 
smooth, corresponding to a roughness spectrum cr(k) 
sharply peaked around k = 0. In this case, we may re- 
place G(k) by G(0) in fgft to find 



SE 



pp 



G(0)(hj 



h\) 



(PFA) 



(30) 



in agreement with Eq. © provided that the response 
function satisfies the limit 



G(k -» 0) 



E'ML) 



(31) 



with the Casimir energy in the ideal case given by 



Epp(L) = HA 



2tt 



d 2 k 



lnfl 



(k,£)V 



-2kL\ 



(32) 

As a consequence of general properties of the rough re- 
flection coefficients at zero momentum transfer (specular 
limit), we show in appendix iBl that this limit is satisfied 
by any response function derived from l|26l) regardless of 
the model considered for the material medium. 



IV. ROUGHNESS RESPONSE FUNCTION FOR 
THE PLASMA MODEL 

In this section, we present an explicit computation of 
the response function G(fc), starting from the general re- 
sult given by Eqs. (|26|I - H29I) . and taking the plasma model 
to describe the optical properties of the metallic mirrors. 
The dielectric function is given by 

LOp 

The plasma wavenumber, wavelength and frequency are 
related by 

kp = — = — . 

Ap c 
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We also define 



Kt (k,0 = \/k 2 +e 



e- 



/bp 



representing the imaginary part of the z component of 
the wavevector inside the metallic medium. The specular 
reflection coefficients are given by 



For given values of k and k', we cast the four coeffi- 
cients App, (£; k, k') into the 2 x 2 matrix A W (k, k') (with 
the association TE = 1, TM = 2), whose nondiagonal el- 
ements represent the coupling between TE and TM po- 
larizations. We find 

A (1) (k, k') = A a) (k, k') - A^ (k, k'), 



r TE (k,0= Kt K 



r™(k,0 = 



K t + K ' 



1 + "|f ) K - Kt 
1 + |l ) < - K, 



(33) 



In order to compute the roughness reflection coeffi- 
cients, we follow the perturbation approach of Ref. [T^| . 
which is based on the extinction theorem [2(J and the 
Rayleigh hypothesis. The incident, reflected and trans- 
mitted fields are related by two integral equations, which 
are solved up to second order of Hj (k) for the reflected 
field in terms of the incident field. This allows us 
to derive the non-specular coefficients R^J, (k, k'; £) and 

Rp (k, k; £) defining the relevant matrix elements of the 
first and second-order reflection operators. 

However, it turns out to be simpler to first calculate 
the coefficients A^l (£; k, k') defined as follows: 



4p'(k,k';£) 



r"(k,Qff'(k',Q (1) 

tP(k,£) V^'W' (35) 



where t p (k, £) are the transmission coefficients for the 
plane interface (see Appendix A). 



A^(k,k') = (Kt±n)B 



-l 
t 



C S 

s q±2M 
' i±ppt i±p/3 t 



(36) 

(34) C = k • k'/(k k') and S = VI - C 2 . We have also defined 



P = - A = - 



B, 



1 
2& 



Primed quantities are likewise defined in terms of k'. 

The second-order coefficients are written in a similar 
way, with 

i?( 2 )(k,k',f) = ^(k,£)A( 2 >(£;k,k'). (37) 
At E (C; k, k') and Atm(£; k, k') are the diagonal elements 
of the matrix 



A^(k,k') = 2kk£- AW(k,k')AL ) (k,k' 



(38) 



with I denoting the 2x2 identity matrix. 

By replacing these results into H28[) and 129|) . we find 



the explicit expressions for the functions anc ^ 



^(0= \ E 



e,e'=+, 



/ TE (k,0 / TE (k',0 C 2 (1 + e/?&) (1 + e'p'ft) + / TE (k,£) / TM (k',£) S 2 (1 + e/3/3 t ) (39) 



+/ TM (k, / TE (k', £) S 2 (1 + ^) + /™(k, £) / TM (k', (C + e0 t ) (C + e'(3'Pt) 



e=4v 



/ TE (k,0 (1-C 2 ^{) (l + ej9A) 



(40) 



+ / TM (k, S 2 (I - /3'#) + / TM (k, (C + e0 t ) (C - ph) 



with 



ft ± K t 

1±/3A 



These expressions can now be applied to the numerical 
computation of the response function for arbitrary values 
of L and Ap. In Fig. 3, we plot G/Epp as a function of 
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0.1 0.2 0.3 0.4 
k [nm" 1 ] 



FIG. 3: Variation of G/Epp versus k for the distances 
L = 50nm (solid line), L = lOOnm (dashed-dotted line), 
L = 200nm (dotted line), and L = 400nm (dashed line). We 
take Ap = 136nm. 



k for several different values of the distance L, and for 
Ap = 136nm, which corresponds to gold covered mirrors. 

According to 125(1 . this ratio provides the relative cor- 
rection of the Casimir energy in the plane-plane config- 
uration when integrated over the roughness spectrum 
<r(k). Moreover, from 10 it also provides the relative 
force correction A for the plane-sphere geometry when 
the sphere radius is sufficiently large. Fig. 3 indicates 
that the relative correction is larger for shorter distances. 

The behavior of G(k) as k — ► is related to the PFA, 
and was already discussed in connection with Eq. ((31(1 . 
In addition to the verification of the general G(k) given 
by ((26(1 - 1(29(1 (see Appendix B), we have also verified inde- 
pendently that the explicit result derived from ((39(1 and 
iflUfl also agrees with Eq. (Jjnj . 

The fact that G increases as k grows from zero, as dis- 
played in Fig. 3, implies that the PFA underestimates 
the roughness correction. In order to quantify the depar- 
ture from the PFA description, we define the sensitivity 
function 



p{k) 



G{k) 

G(oy 



(41) 



In Fig. 4, we plot p as function of k for the same values of 
distance and Ap employed in Fig. 3. The PFA amounts 
to replace p(k) by unity for all values of k contained in 
the roughness spectrum a(k). Clearly, this approxima- 
tion is better for shorter distances, and smaller values 
of k (corresponding to longer roughness wavelengths) , as 
expected. For instance, the inlet shows that the PFA is 
a good approximation for L = 50nm and k < 0.04nm _1 . 
On the other hand, for L — 200nm and k — 0.02nm 



(roughness wavelength 2ir/k ~ 300nm) we find p ~ 1.6, 
corresponding to a roughness correction 60% larger than 
the PFA result. 




0.0 0.1 0.2 0.3 0.4 
k [nm 1 ] 

FIG. 4: Variation of p versus k for several values of L (same 
conventions as on Fig. 3). 



Fig. 4 also indicates that p(k) grows linearly for large 
values of k. In the next section, we show that this is a 
general result, valid for arbitrary values of L and Ap. 



V. HIGH-fc LIMIT 

When the momentum transfer k is much larger than 
1/L, the function ^ k '_ k (£)j representing the contribu- 
tion of the first-order reflection operator, is negligible, be- 
cause it is proportional to the exponentially small prop- 
agation factor exp(— ^/(k' — k) 2 + t; 2 L) w exp(— kL) ap- 
pearing in Eq. ((28(1 . This general property can be un- 
derstood from the discussion of Sec III: the loop function 
<5/p 2l ^(k', £) contains two rough reflections separated by 
a intracavity round-trip propagation with the modified 
momentum k' — k. 

On the other hand, the loop function Sfp 2u \li, £) in- 
volves a single second-order rough reflection, which must 
conserve momentum so as to allow for a closed loop. 



Thus, b 



(ii) 
k',k' 



does not involve propagation with the 



modified momentum k' — k and is the dominant term 



in the high-k limit. We calculate b 



(ii) 
k'.k' 



from (gOl by 



taking k' — k w — k. We also assume that k 3> fcp; the 
opposite case will be discussed in Sec. VI. We find 



p(k)~ak for k 1 ^;Ap,L. 



(42) 



The dimensionless parameter a/L depends on Kp 
kpL = 2nL/Xp only, and is given by 
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hcA r°° f, Kl 



K 2 



7/ TE (k,0 



2( 7 2 -fi 2 ) 2 - 7t 2 (2 7 



3ft 2 



{lit) 2 - (l 2 ~ ft 2 ) 2 



-7/ TM (k,£) 



(43) 



We have introduced the dimensionless integration vari- 
ables 7 = kL, -ft — K t L, and ft = £L/c. / TE (k, £) and 
/ TM (k, £) ar e calculated from lfTT)|) as functions of 7 and 

ft. 

We plot the coefficient a as a function of L in Fig. 5, 
with the plasma wavelength of gold Ap = 136nm as in 
the previous numerical examples. At the limit of short 
distances, we recover from (I43|l our previous result |Tlj 



a = 0.4492L for k^ 1 <L< A r 



(44) 



This corresponds to the high-A; limit of the plasmon 
(non-retarded) regime, which we shall discuss further in 
Sec. VII. This limit is indicated by the dotted line in 
Fig. 5. 

As shown by Fig. 5, the angular coefficient a saturates 
at the limit of large distances. This corresponds to the 
limit fc _1 <C Ap <C L, which may be obtained analytically 
from Eq. (|43l) by expanding its r.-h.-s. in powers of Ap. 
The integrand vanishes to order Ap 1 , whereas the the 
zero-th order term yields 



hcA 



d-f- 



1 



(2tt) 2 L4G(0) Jo 

[ dft( 7 2 - 3ft 2 ) = 
Jo 



27 



(45) 



so that the dominant term is of the order of Ap. We also 
need to calculate G(0) in the limit L ^> Ap. As expected, 
we find G(0) = E$ T p "(L)/2 = - tt 2 hcAj (120 L 5 ), where 
Epp is the Casimir energy for perfectly reflecting mirrors. 
We then find 



60 X f A 
a =— Ap / ar/ 

Jo 

P 

x 



7T 



=2 7 



7 ( e 27 - 1)2 

f dn(~f 4 ~ 2j 2 n 2 + 3ft 4 ) 

Jo 



(46) 



giving 



15tt 



-Ap for 



< A P < L, 



(47) 



which is in agreement with the saturation value shown in 
Fig. 5. This result remarkably differs from the long dis- 
tance behavior reported in Ref. [TT| | . which corresponds 
to the perfectly- reflecting limit. Note that the high- A: ex- 
pression (|43Jl holds when the roughness length scale 1/fc 
is much smaller than both Ap and L. In this regime, and 
as a consequence of the momentum transfer induced by 
the roughness effect, the modified field momentum has a 
magnitude |k' — k| much larger than fcp. Therefore, it is 
poorly reflected by the mirrors, even though the initial 




200 400 600 800 



L [nm] 



FIG. 5: Variation of the angular coefficient a versus L for 
Ap = 136nm. The analytical result for fc" 1 <C L C Ap is 
shown as the dotted line and for fc _1 C Ap < I as the 
dashed line. A comparison between this second result (dashed 
straight line) and the exact p(k) (solid line) is shown in the 
inlet for L = 2/im. The analytical result p = Lk/3 predicted 
by the model of perfect reflectors (dotted line) is valid only 
in the intermediate range Ap < < L. 



momentum satisfies k' ~ 1/L <C kp. In order to obtain 
the perfectly-reflecting limit, one must assume that Ap 
rather than 1/k is the shortest length scale, as discussed 
in the next section. 



VI. PERFECTLY-REFLECTING LIMIT 

In this section we assume that Ap is much smaller than 
both the separation L and the roughness wavelength 1/k. 

rep- 



(ii) 
k',k" 



Then, we expand 6k'. k" in powers of Ap. For b 
resenting the contribution of the second-order reflection 
operator and given by (|40|l . the dominant term is of the 
order of 1/Ap : 



6 W _ 47 r«lk' 
k K k k' Ap 



k + O(\ p) , 



(48) 



where / pr (k') represents the loop function for perfect 
reflectors [it is the same for both polarizations, with 
r TE = _ r TM = _ lj accor ding to @ and l|3ljl]. On 

the other hand, b^'V-k vaiusnes up to order 1/Ap. ft 
follows that G (k) vanishes at this order, because when 
taking the integral of the r.-h.-s. of (|48|l over all values of 
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momentum k' in l|26|l , opposite values of k' compensate 
each other. 

Both fr^k" and k" contribute up to order A P . It is 
useful to replace 6k', k" by the symmetrized form 

bk'M" = (&k', k" + 6 k ",k')/2. 

This procedure does not change the response function 
because G(k) = G(-k) = G(k). Taking d 2 k' = dk'k'dcf)' 
in 1)26(1 . we derive 



dk' k' 



b' b 



k'.k' 



(49) 



6k', k' 
gDJ: 



is obtained from the term of order A P in 1(39(1 and 



-2k' L 



-2k" L 



6k', k" = 



(1-e 



-2k' L 



)(1 







(50) 



(k'k") 2 + (£ 2 /c 2 + k' -k")' 



We change the variables of integration from (£, fc') to 
(ft', k"). The integral over <// yields 



|J| k{K',K",<j>')d4>' 



7T ft ft 

2~F 



where </> m = arcsin[(ft' 2 — k" 2 + fc 2 )/(2ft'fc)] and |J 
Jacobian corresponding to the transformation. 
From l(l^|) - (|^T|l we derive 



(51) 
is the 



G(k) 



1 



hcA 

~ 8^ Uq J 1 

(77') 2 + |(7 2 +Y 2 



d'ye 



-27 



e -2 7 
„2\2 



7+9 



di (52) 



1 



-2 7 ' 



It— si 
for Ap 



0. 



7 has the same meaning already discussed in connection 
with l(43|) while 7' corresponds to the diffracted wave. 
From O we verify that G(0) = E^ p "/2 as expected. 
For arbitrary values of q = kL, numerical i nteg ration 
of 152() agrees with the results of Emig et al. 21] for a 
perfectly-reflecting mirror corrugated along a fixed direc- 
tion in the xy plane. By taking the assumption of perfect 
reflectivity from the start, we may derive this result for a 
general deformation directly from our general expressions 
(|2*6*)l - lj2l?)l . as discussed in Appendix C. 

In order to discuss the regime Ap <C 1/fc <C L, we 
now take the high-fc limit of the right-hand side of (152(1 . 
Due to the presence of the exponential factor exp(— 27), 
the dominant contribution comes from the corner 7 ~ 
1, 7' ~ q of the rectangle associated to the integration 
region. We may thus neg lect exp(— 27') and recover the 
long distance limit of 



G{k) = -- 



2 HcAq 



7 3 g -2 7 



e- 2 7 360 L 4 ' 



-Lk for A P < k 1 < L. 



(53) 




10 2 10? 
V[ nm ] 



FIG. 6: Variation of p versus Ap for k = 0.02nm and the 
distances L = 2/im (solid line), L — 600nm (dashed line) and 
L — lOOnm (inlet, solid line). The values at Ap = 136nm are 
indicated by vertical lines. 



In summary, the long-distance behavior is given by 
(IT1 when 1/fc < A P < L, and by (JHSl when A P < 
1/fc L. The cross-over between these two regimes 
is shown in the inlet of Fig. 5, where we plot p as a 
function of fc for L = 2/im. The finite conductivity of 
the metals clearly reduces the roughness correction for 
very large values of kL, due to the saturation effect dis- 
cussed in Sec. V. In Fig. 6, we plot p as function of 
the plasma wavelength Ap for k = 0.02nm _1 (rough- 
ness wavelength 27r/fc « 300nm), with L — 2pm (solid 
line) and L = 600nm (dashed line), in order to ana- 
lyze in detail the effect of finite conductivity. The value 
Ap = 136nm corresponding to gold covered mirrors is 
highlighted by a vertical line. The perfeclty-reflecting 
limit corresponds to the small values of Ap shown in the 
left-hand side of the figure. The correction decreases as 
Ap approaches 27r/fc because of the saturation effect, and 
then it increases again as Ap approaches and goes be- 
yond the separation distance L. For Ap 3> L, we find the 
limit predicted by the plasmon model, to be discussed 
in the next section, which is larger than the perfeclty- 
reflecting limit llj. When kL ~ 1, the correction is 
always larger than the perfectly-reflecting limit, as illus- 
trated in the inlet of Fig. 6, where we take fc = 0.02nm _1 
and L = lOOnm. In this case, the perfectly-reflecting 
limit is 11% larger than the PFA result, whereas at 
Ap = 136nm we find a correction 19% larger than the 
PFA. 
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VII. PLASMON LIMIT 



In the short distance regime, L <C 
energy is associated to surface plasmons 
the previous sections, when integrating 
the dominant contributions come from 
£ such that n' ~ 1/L. Then, the short- 
ly yields r TE = 0(L/A P ) 2 , so that 
of TE polarization is negligible. From 
negligible except for the values 

£ ~ wp -C c/L, 
for which k rs k , and 



Ap, the Casimir 
. As discussed in 

&k',k'-k in ©, 
values of k' and 
distance limit of 
the contribution 
(EH, r™ is also 



2£* + w* 



(54) 



(55) 



In other words, the dominant contribution is associ- 
ated to low-frequency surface waves, the reflection co- 
efficient having poles at the surface plasmon resonance 
£ = -wp/y/2. From we also conclude that the re- 
tardation time L/c for propagation between the mirrors 
is negligible in the time scale associated to the relevant 
field frequencies. 



We calculate the roughness correction by taking the 
appropriate limits in ((39(1 and 1)40(1. Since r™ only de- 
pends on £, the diffracted wave sees the same reflection 
coefficient, which we denote as r for simplicity. We find 



fP - 

k '' k " " 2(1 



-2k' 'L g-2fc" ' L 



r 2 e- 2k ' L ){\ 



j,2 g — 2/c" L \ 



[(C + l) 2 + 2r(l - C 2 ) + r 2 (l - C) 2 



(56) 



°k',k" 



-2k' L 



re 



-2k' L 



[r(l - C) 2 + 1 - C 2 



(57) 



These equations do not agree with the results of Ref . ^(| , 
which in their turn are not consistent with the PFA. 
Therefore, it is important to check our results against 
the PFA by taking k" = k' = k and C = 1 in and 



&k,l 



2A .2 r 2 e -2fci 



(1 



(58) 



When replacing this result into (|26|) . and taking l|55|) 
into account (with £ = £/k>p, and 7 = kL = fc£), we 
obtain 

C(0) = fdcft, 2^g±ig , 

2^ 2 L 4 y Jo [(2C 2 + l) 2 e 2 ^ - l] 2 

Integrating this expression by parts, one shows that G(0) 
satisfies 1(31(1 . with Epp representing the short-distance 
limit of the Casimir energy as given by : 



Epp(L) 



HA uip 

47T 2 L 2 



00 7 2 

dC /o d7 (2C 2 + l) 2 e 2 ,-F 



This confirms consistency with the PFA, which is in line 
with the more general discussion presented in Appendix 
B. 

In addition to the limit k — > 0, it is also interesting to 
analyze the case k ^> 1/L from ((56(1 and ((57(1 . allowing for 
the evaluation of the limit 1/i < L < Ap. This provide 
as additional check, by comparison with the results of 



Sec. V. In this limit, b k / k ,_ k is exponentially small, and 

the dominant contribution for 6 k / k >_ k comes form the 
second term in (|57|) , which is proportional to k : 



^k',k'-k ~ &k', -k 



! ; e -m [ki - c) 2 + 1 - c 2 ] *. 

(59) 

We also take C m — k' • k/(fcfe') = — cos</>, and integrate 
over </> to derive from 1(26(1 



(60) 



7 2 (2C 2 +4) 



X / dC X d7 (2C 2 + l) 2 [(2C 2 + l) 2 e^-l]' 

yielding p(k) = 0.4492L k, in agreement with Eq. l|4^| 
and Ref. Fnl. 



VIII. CONCLUDING REMARKS 

We have calculated the second-order response func- 
tion G(k) for arbitrary values of the plasma wavelength 
Ap and the distance L. This allows for a reliable com- 
putation of the roughness correction, once the roughness 
spectrum cr(k) characterizing the metallic surfaces is ex- 
perimentally determined. In order to gain further insight 
into the the roughness correction itself, we consider the 
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particularly simple example of a Gaussian spectrum ^(| : 



cr[k] = 7ra 2 ^exp 



k 2 ^ 



(61) 



10" 



The roughness variance values a 2 , and £c represents the 
correlation length. 

According to Eq. I|25|l . the relative force correction 
A is obtained by integrating the normalized response 
function G(k)/Epp (see Fig. 3 for some numerical ex- 
amples) over the Gaussian spectrum given by (|61|l . In 
Ref. [l4[, we have discussed some simple analytical ex- 
pressions in the limiting cases Xp <§; tc and Xp 3> £c, 
which can be easily derived from the results of Sees. V- 
VII. However, the experimental parameters are likely to 
be such that neither of the two limits holds. Hence, we 
must rely on the numerical calculation of G{k) to com- 
pute the correction. In Fig. 7, we compare the exact 
results for A/a 2 (for two different values of £c) with the 
PFA formula A/a 2 = E^ F (L)/[2E PP (L)\, taking as be- 
fore Xp = 136nm. At L = lOOnm, the exact value is 57% 
and 7% larger than the PFA result for £c — 50nm and 
lc = 150nm, respectively. In agreement with the discus- 
sion of Sec. IV, the PFA is better for shorter distances 
and longer correlation lengths. For instance, as suggested 
by Fig. 7, it provides accurate results if L < lOOnm and 
£c > 150nm. Note, however, that the validity of the PFA 
can be addressed in a reliable way only from the analysis 
of the experimentally measured roughness spectrum. In 
the Gaussian model discussed in this section, the contri- 
bution of high values of k are exponentially small, but in 
the real case the decay of a(k) might be smoother. Thus, 
this approximation might be worse than discussed here. 

For large values of L, no simple analytical result is 
available when Xp ~ lc- In the inlet of Fig. 7, we show 
that the perfectly-reflecting result A/a 2 = l^pn / (£cL) 
valid for Xp <C £c <C L overestimates the correction for 
£c = 50nm by more than 50%. This is a consequence of 
the saturation effect discussed in Sec. V: diffraction by 
roughness Fourier components at k > Xp 1 give rise to 
waves which are poorly reflected by the mirrors. Thus, 
when the roughness spectrum contains very high values 
of k, the correction is reduced with respect to the result 
of the perfectly-reflecting model. The reduction factor 
decreases as kXp — ► oo down to the limit (7/5tt)Xp/L <C 
1. 

In conclusion, we presented a perturbative method for 
the calculation of the Casimir energy between rough mir- 
rors, up to second order in the amplitude of the deforma- 
tion. It relies on the manipulation of reflection operators, 
taking into account diffraction and the coupling between 
different field polarizations. We applied the method to 
compute the roughness correction in the framework of 
the plasma model, for arbitrary values of the plasma and 
roughness wavelengths and the mirror separation L. An- 
alytical results for different limiting cases were discussed. 
In particular, the PFA regime follows from our formalism 
in the limit of very smooth surface profiles. By compar- 
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FIG. 7: Roughness correction for the Gaussian spectrum. We 
plot the relative force correction for the plane-sphere setup 
over the squared amplitude of roughness, A/a 2 , versus L for 
lc = 50nm (dashed line) and la — 150nm (dotted line). The 
solid line represents the PFA result. In the inlet, we also plot 
the long-distance perfectly-reflecting limit (dot-dash). We 
take Ap = 136nm. 



ison with our numerical results, we were able to analyze 
the accuracy of the PFA in the problem of roughness. 
For a given roughness spectrum, our theory provides re- 
liable numerical results for the roughness correction, and 
allows to check the validity of the PFA approach in a 
given experiment. More realistic models for the metallic 
mirrors |7| can also be considered by applying the formal 
results presented here. 

We thank Cyriaque Genet and Marc-Thierry Jaekel 
for discussions. PAMN thanks Instituto do Milenio de 
Informacao Quantica and CNPq for partial financial sup- 
port. 



APPENDIX A: OPTICAL NETWORK THEORY 
AND THE CASIMIR FORCE 



In this appendix, we compute the spectral density in 
the intracavity re gion by generalizing the optical network 
formalism of Ref. [13| to the case of rough surfaces. This 
allows us to compute the Casimir force with the help of 
the Maxwell stress tensor. 

The basic idea is to derive the relation between the 
intracavity field and the incoming outside field, whose 
fluctuations are known. The outside field propagating 
from the region z < (see Fig. 1) is written as [with 
r = (x, y), u> = cy/ k 2 + k 2 z and eo denoting the vacuum 
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a column vector: 



FIG. 8: Input and output fields. 



permittivity] 



p-TE,TM 



(Al) 



x exp[z(k • r + k z z — ujt)]e p + H.c. 



The Fourier components e™ p (K) satisfy the commuta- 
tion relations of freely-propagating fields: 



;K),er (K') 



(27r) 3 5 (2) (k-k')(5(fc 2 



_(A2) 

The free-space fields propagating along the negative z- 
direction are written in a similar way, except for the re- 
placement exp(?fc z z) — > exp(—ik z z). We use the same 
K — k + k z z to label them as well, and our notation is 
such that k z > in all cases. 

We define scattering and transfer operators for the two 
rough mirrors and for the empty-space propagation be- 
tween them. The cavity is taken as a composed net- 
work, and the corresponding transfer operator is simply 
the product of the transfer operators for the elementary 
components. 

The scattering operator characterizing a given element 
of the network is defined in the following way. As shown 
in Fig. 8, the input field contains components propagat- 
ing from the left and righthanded sides, and with polar- 
izations TE and TM. We arrange these components into 



e in (K) 



„inTE 

ilITM 
„hlTM 

e R 



/k 



We employ a similar notation for the output field: 



The scattering operator provides the input-output re- 
lation 




fc (K) = 



d 3 K' 



S(K,K')e in (K') 



(A3) 



i( K >0) (2tt) 3 
where the integral over k' runs from to oo. 

Whereas Ref. 13] allows for lossy mirrors, here we as- 
sume that there is no dissipation in our network. Hence 
we only consider unitary scattering operators: S • S^ = 1, 
where 1 is the identity operator. In terms of the corre- 
sponding matrix multiplication, this condition reads (I is 
the 4x4 identity matrix) 



d 3 K' 



(27T 



gSCKi.K'JSCKa.K'jt = (27r) 3 <^ 3 )(Ki— K 2 ) I. 

(A4) 

We analyze in detail the scattering operator Si corre- 
sponding to mirror Ml in Fig. 1. We expand Si up to 
second order of the deformation amplitudes of the two 
lateral mirror's rough surfaces: 



Si 



j(0) 



SS 



(2) 



(A5) 



The scattering by ideal plane surfaces conserves to, k and 
polarization: 



S { °> (K, K') = (2tt) 3 5 (3) (K - K') 



( *T"(k,w) 
^»(k,w) 




4-TE 

n 



\ 







(k,w) 
(k,w) 





J.TM 

E l 
' 1 






(k,w) 
(k,w) 



,.TM 
1 



\ 



(k,«) 



r^(k, u>) are the specular reflection coefficients as seen by the intracavity field presented in Sec. III. The coefficients 
^(k, oj) play the same role for mirror M2. Together they turn out to be the only relevant ones for the calculation 
of the Casimir effect, both in the ideal and rough cases. For simplicity, we have denoted rf(k, u>) simply as r^(k, u>) 
everywhere in this paper, except in the present appendix. 

The first (£ — 1) and second-order (i = 2) corrections mix up different polarizations and values of k, but conserve 
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the frequency: 



SS[ e) (K,K') = 2nS(k' z - \A 2 - k' 2 + k 2 



( 


TE;TE] 


6r® 


[TE;TE] 


5%> 




[TE;TE] 


5tf 


TE;TE] 


Sr[ e) 


sap 


TM;TE] 


8r<P 


TM;Te] 


5^ 




TM;Te] 


8tf 


TM;TE] 





[TE;TMJ 
[TE;TM] 
TM;Tm] 



TE;TM 



Stf^ [te ; tm] 
5rf^ [tm ; tm] 



(A6) 



All matrix elements above are functions of frequency 
and of the initial and final momenta k' and k. 
The four elements 8rf^ [p; p'] are the matrix elements 

(k,p| 572.^ \k',p') introduced in Sec. III. 

The transfer operators T provide the fields at the left- 
hand side of the mirror in terms of the fields at the right- 
hand side. They can be obtained from the scattering 
operators as follows [l3| : 



t = - (p_ - s p+r 1 • (p+ - s • p_) . 

The operators P+ and P are defined by 



P+(K, K') = (27r)^ J (K - K') 

P = 1 P + . 

The intracavity field (see Fig. 1) 



/ 1 



10 

\ 



e°(K) 



G 

TM 



\ e c J 



K 



(A7) 



is computed from 



e°(K) 



(*;>o) 



d 3 K' 
(2^)3 



R(K, K')e m (K') 



with 



P — T2 • P+ • S cav -(- T2 • P_ , 



(A8) 



(A9) 



and where S ca v is the scattering matrix for the cavity as 
a composed network. 

The Casimir force on mirror Ml is computed from the 
energy-momentum tensor component T zz , evaluated at 
the intracavity region and at the outer side of the mirror, 
and averaged over the vacuum state: 



^pp = / d 2 r [<T 2 L z (r)) vac - (T£ (r)) vac ] . (A10) 



We calculate T^ z (r) using Eq. i|Al|l and similar expan- 
sions for the magnetic field and for the outgoing 
fields: 



1 f d 3 K r°° dk' 

d 2 r<r z L z (r)) vac = - / —-3 / ^^cos 2 ^([er(k,A: z )er(k,^)t + er p (k,fe)er p (k,fc;) t ]), 

2 J(k,>o) ( 27r ) Jo 27T _ 



>(k z >0) 



P 

, (All} 



with cos 8 = k z /K. A similar expression is found for the 
inner region in terms of the intracavity fields. 

When taking the average over the vacuum state we use 
the commutation relation l|A2(l to find 

(e£*(k, k z )er' (k, KY) vac = 2n A S(k z -k' z ). (A12) 

The outgoing field e£ ut satisfies the same commuta- 
tion relation, and provides an identical contribution in 
Eq. (XH) . 

On the other hand, the commutation relation of the 



intracavity field is modified by the joint effect of the two 
mirrors. We derive the corresponding spectral density 
from Eqs. 0, (S]^ and (jA"l2|) : 

(ec (k,k z ) 7 P C (k,^)t) vac = 2nAg p {k,uj)8(k z - k' z ), 

(A13) 

We obtain the explicit expressions for Sfj, 21 ^ (k, lo) and 

6/^ 2i) (k,w) given by Eqs. O and $M with the hel P 
of a computer algebra system. Remarkably, they only 
contain the rough reflection coefficients associated to in- 
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ternal reflections [only 4 out of 16 elements in Eq 
for instance]. 

The spectral density for e q is modified, with respect to 
the free-space case, by the same generalized Airy function 
<7p(k, u>), as far as second order terms containing rough 
reflections at the same mirror are concerned. After re- 
placing HA12jl into IjAlip . and using (|A10|) and (|A13j) . we 
derive the result given by Eq. I|14(l for the Casimir force. 



APPENDIX B: PROXIMITY FORCE 
APPROXIMATION AS A LIMITING CASE 

In this appendix, we derive the PEA result for the en- 
ergy correction as a limiting case of the general results 
of Sec. III. This will bring a deeper understanding of the 
PEA, by showing its connection with some specific prop- 
erties of the reflection operators. These properties must 
be satisfied regardless of the particular model considered 
for the material medium, and are related to the specular 
reflection by displaced plane mirrors (specular limit). 

In the PFA regime, the Fourier profile functions 
Hj(Ak) are sharply peaked around Ak = 0. Thus, we 
may replace k' by k in the argument of the non-specular 
coefficients appearing in the r.-h.-s. of l|2UI) and (|21|l . The 
resulting expressions correspond to the case of ideal plane 
mirrors which are displaced from the associated reference 
planes at z = and z = L. In this case, the reflection is 
modified with respect to the non-perturbed case just by 
the effect of the propagation from the reference plane to 
the mirror and back. Up to second order, this amounts 
to take the series expansion of the exponential factor 



yielding 



and 



^(k,k;0 = -2«rJ(k,0« PJ 



flg(k,k;0 = 2 K 2 7j(k,0. 



(Bl) 



(B2) 



Eqs. (|B1I) and l|B2|) are general properties of the reflection 
coefficients, and are useful for checking explicit calcula- 
tions, regardless of the specific model considered for the 
material medium. When using these results to compute 
G(0) from Eqs. H>|)-(H5J), we obtain Eq. (J2H, with the 
Casimir energy in the ideal case given by l|32[l . This ver- 
ifies the PFA limit as discussed in the end of Sec. III. 



APPENDIX C: PERFECT MIRRORS 

In Sec. V, the perfectly-reflecting limit was derived 
from the plasma model results by taking Ap <C k~ l ,L. 



This appendix presents an alternative, simpler deriva- 
tion, in which the usual model of perfect reflectors is 
taken from the start. The case of corrugation along a 
fixed direction in the xy plane (say the direction along 
the a;-axis, with h\{x,y) — h\(x)) was considered by 
Ref. [2l|. In this case, the calculation can be consider- 
ably simplified by taking a convenient definition for the 
field polarizations, which then turn out to be not coupled 
by the scattering from the surface. On the other hand, 
in this paper we consider arbitrary small-amplitude de- 
formations, so that the coupling between different polar- 
izations has to be taken into account. 

For the mirror near z — 0, we take the boundary con- 
dition 



ni(a;, y) x E(x, y, hx{x, y)) = 0, 



(CI) 



where h\ (x, y) is the unitary vector normal to the tangent 
plane at the point {x, y) : 



ni(x,y) 



V/zi 



v/i + (vM 2 



As explained in Sec. II, positive values of hi(x,y) are 
defined along the positive z axis, and the (intracavity) 
incident field, as given by J7J propagates along the neg- 
ative z axis with a phase factor e~ tkzZ . 

Solving I|C1|I up to first order of hi(x,y), we deter- 
mine the complete first-order reflection operator. Using 
the matrix notation introduced in Sec. IV, the first-order 
non-specular coefficients defined by Eq. (|21|l are written 
as 



R (1) (k,k';w) = 2i 




(C2) 



where C and S are the cosine and sine of the angle be- 
tween k and k'. 

We only need the diagonal elements of the second-order 
reflection operator, which we collect from the solution of 
(jClf) up to second order. The corresponding non-specular 
coefficients [see Eq. (|2(J|) ] are 



4 2 2(k,k'; u) = 2k z k' z C 2 + 2^S 

c z k' 



(C3) 



MS(k,k'; W ) 



k,k'. 



(kk' -uj 2 C/c 2 Y 



^ 2 K n2 



(C4) 

These results satisfy the specular limit discussed in Ap- 
pendix B. When replaced into the general expressions 
of Sec. Ill, they reproduce the results for the perfectly- 
reflecting limit of Sec. V. 
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